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Abstract. Rigourous calculations of the imaging properties of metamaterials 
consisting of metal-coated semiconductor nanoparticles are presented. In particular, it 
is shown that under proper choice of geometric and materials parameters, arrays of such 
particles exhibit negative refractive index within the region of the excitonic resonance 
of the semiconductor. The occurrence of negative refractive index is predicted by the 
extended Maxwell-Garnett theory and confirmed by a layer- multiple scattering method 
for electromagnetic waves. By using the same method it is shown that within the 
negative refractive-index band, arrays of such nanoparticles amplify the transmitted 
near-field emitted while simultaneously narrow down its spatial profile leading to 
subwavelength resolution. The effect of material losses to the imaging properties of 
the arrays is also addressed. 
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1. Introduction 

According to a recent definition [1], metamaterials are artificial materials which 
exhibit response characteristics that are not observed in the individual responses 
of its constituent materials. The most fascinating class of metamaterials are those 
exhibiting simultaneously negative permittivity e and permeability /i, i.e., a negative 
refractive index (NRI) [21 E] • The usual approach of designing a NRI metamaterial is 
to combine a 'magnetic' sublattice (one which exhibits negative /i) of miniaturized RLC 
circuits, e.g. split-ring resonators, with an 'electric' one (exhibiting negative e) of thin 
metallic wires [1]. However, an alternative route has been recently suggested where 
the electric and magnetic sublattices are occupied by units of less elaborate geometry 
such as cylinders or spheres made from resonant materials (ionic, semiconducting or 
plasmonic materials) displaying high-refractive index within a specific frequency window 
ilElIZlElElIinillllinillSlEllISlIIS]. The magnetic activity of the cylinders / spheres hes 
within the region of Mie resonances resulting from the enhancement of the displacement 
current inside each sphere which, in turn, gives rise to a macroscopic magnetisation of 
the whole structure. The electric activity is attributed to the large polarization induced 
to the sphere due to the giant dielectric permittivity around the resonance frequency 
of a given dielectric function, e.g. polariton, plasmon or exciton resonance. The above 
theoretical suggestions have been recently verified experimentally for arrays of dielectric 
particles of millimeter [T7] and micrometer [TH] size. Such arrays can be miniaturized 
to the nanometre scale in the form of nanoparticle superlattices [19] . 

Perhaps, the most paramount application of NRI metamaterials is their imaging 
property; as has been pioneered by Pendry [3] and verified experimentally, a planar slab 
of metamaterial with NRI can overcome the standard diffraction limit in imaging by 
focusing the far-field via negative refraction [SUl EI] and by amplifying the near-field 
[221 [23l [2l] by surface-plasmon excitation, a possibility which promises the realization 
of a perfect lens. Near-field amplification is also feasible separately for S-polarized waves 
in the magnetostatic limit using solely RLC circuits such as split-ring resonators and 
Swiss rolls [251 126] and for P-polarized waves in the electrostatic limit using periodic 
structures with metallic components such as rods ^7\, wires [2B], and spheres [29 j . 

In this work we study the imaging properties of three-dimensional (3D) NRI 
metamaterials consisting of metal-coated semiconductor spheres. As has been recently 
shown, arrays of semiconductor spheres with strong excitonic oscillation strength (CuCl, 
CU2O) can exhibit negative permeability in the visible region, around the exciton 
resonance [H]. In order to realize an NRI structure in the regime of negative /i, the 
array of semiconductor spheres can be (a) either combined with an array of metallic 
spheres (negative e) in a way that a binary NRI structure is created [8j, or (b) by 
coating the semiconductor spheres with a metallic shell [IIl|30]. The above structures 
have been studied by employing both effective-medium treatments [HI [121 [13] and more 
rigorous electromagnetic approaches [H [H] . By exploiting the occurrence of NRI in the 
case of arrays of metal-coated semiconductor nanoparticles, we will show that near-field 
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amplification along with suhwavelength resolution of an image can be achieved in such 
structures, for both polarization modes. We will also elucidate the role of absorption 
which is intrinsic in the materials constituting the NRI metamaterial under study. The 
paper is organized as follows. In section [2] the basic theoretical tools for studying the 
imaging properties of arrays of spherical inclusions is presented. Section [3] applies the 
theory to the case of metal-coated semiconductor spheres and section H] concludes the 
paper. 



2. Theory 

In this work, we are dealing with finite slabs of metamaterials consisting of a number of 
planes of spheres with the same 2D periodicity. In order to probe the imaging properties 
of these structures, we consider the electric field emitted by a localized source, namely, 
that radiated by a small, linear, infinitely thin, center-fed antenna. The antenna lies 
along the 2;-axis from —d/2 < z < d/2 with current /. We expand the electric field 
radiated by the antenna as a series of spherical waves [31] 

oo I f . 

= HY^ \amraht{qr)y.i^{i) + aEim-V x [/i+(gr)Xi^(f )] i , (1) 

1=1 m=-l ^ 1 ) 

Xim(f) are the so-called vector spherical harmonics j31] and hi are the spherical Hankel 
functions of order I. q = cu/c, where c = 1/y/JIeJIoeo = cq/ y/JIe is the velocity of light 
in the medium surrounding the antenna. The magnetic multipole coefficients anim are 
zero whilst the electric ones are given by [32] 

_ / l 4n{2l + l) qd , qd 

provided that / is odd and m = 0. jz in equation denote the spherical Bessel 
functions. Note that the above formula is valid in the case where the antenna is much 
smaller than the wavelength [32], i.e., qd <^1. 

Since we wish to study the transmission of the above field through a slab of a number 
of periodic planes of spheres, it is advantageous to transform the field of equation ([T]) 
to a basis of plane waves consistent with the 2D periodicity of the planes of spheres. If 
the linear antenna is placed to the left of the slab (see the calculation setup in figure H]), 
then the field radiated to the right and is incident on the slab is written as (assuming 
the center of coordinates be located at the localized source) [33] 

+(r) = ^ / / d\ Yl Er exp(iK+ ■ r) (3) 

^0 J J SBZ g 

with 

oo I 

^(k||) = E E E ^Pim-AKypim (4) 

1=1 m=-l P=E,H 
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where i = 1,2 are the two independent polarizations (polar and azimuthal) which are 
normal to the wavevector [311 ESI ES] 

i/r 



K: 



-(k|l+g)' 



(5) 



The vectors g denote the reciprocal-lattice vectors corresponding to the 2D periodic 
lattice of the plane of spheres and ky is the reduced wavevector which lies within 
the Surface Brillouin Zone (SBZ) associated with the reciprocal lattice [SI EHl 136] . 
When < (ky + g)^, the wavevector of equation defines an evanescent wave. The 
coefficients Apim are given by 
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where 6, denote the angular variables (Kg ) of Kg and Aq is the area of the unit cell 
of the 2D lattice occupied by the spheres, l^*" denotes a spherical harmonic as usual, 
ei, 62 are the polar and azimuthal unit vectors, respectively, which are perpendicular to 
Kg. a™ are given by a™ = ^[{1 — 'm){l + m + 1)]^''^. The incident field of equation ([3]) 
will be partly transmitted through the slab under study. The transmitted field will be 
given by 



E 



tr 



^ [ [ rf^A:||^E*'-+(k||)exp[iK^ 

^0 J J SBZ „ 



d)] 



with 



(7) 



(8) 



d is a vector joining the source to the image (see figured]). The transmission matrix 
appearing in equation ([H]) is calculated within the framework of the layer-multiple- 
scattering method which is an efficient computational method for the study of the EM 
response of three-dimensional photonic structures consisting of nonoverlapping spheres 
[Ml [35l [36] and axisymmetric non-spherical particles [37] . The layer-multiple-scattering 
method is ideally suited for the calculation of the transmission, refiection and absorption 
coefficients of an electromagnetic (EM) wave incident on a composite slab consisting of 
a number of layers which can be either planes of non-overlapping particles with the 
same 2D periodicity or homogeneous plates. For each plane of particles, the method 



Suhwavelength imaging of light by arrays of metal- coated semiconductor nanoparticles. . .5 



calculates the full multipole expansion of the total multiply scattered wave field and 
deduces the corresponding transmission and reflection matrices in the plane-wave basis. 
The transmission and reflection matrices of the composite slab are evaluated from those 
of the constituent layers. By imposing periodic boundary conditions one can also obtain 
the (complex) frequency band structure of an infinite periodic crystal. The method 
applies equally well to non-absorbing systems and to absorbing ones. Its chief advantage 
over the other existing numerical methods lies in its efficient and reliable treatment 
of systems containing strongly dispersive materials such as Drude-like and polaritonic 
materials. 

The calculation of the incident [equation ([3])] as well as the transmitted field 
[equation ([7])], requires a numerical integration over the entire SBZ. In the example 
examined in the next section, the spheres in all planes occupy the sites of a square 
lattice and, therefore, the SBZ is also a square. The SBZ integration of equations ([3]) 
and ([7j) is performed by subdividing progressively the SBZ into smaller and smaller 
squares, within which a nine-point integration formula [3_8J is very efficient. Using this 
formula we managed excellent convergence with a total of 73728 points in the SBZ. Also, 
the inclusion of 13 reciprocal-lattice g-vectors along with an angular-momentum cutoff 
imax = 4 provided converged results. 

3. Results 

We consider a 3D array of closed-packed CuCl nanoparticles of radius S = 28 nm; 
CuCl exhibits a Z3 exciton line at 386.93 nm |39]. Around the exciton frequencies, the 
dielectric function of the above semiconductors is given by 

tsi^^) = ^00 + A'y/{uJo-uj -ij). (9) 

The constant A is proportional to the exciton oscillator strength and for CuCl, 
A = 632. The rest of the parameters for CuCl are [32]: = 5.59, Jtluq = 3.363 eV, 
and ^7 = 5- 10~^ eV. The small value of the loss factor 7 implies a very narrow 
exciton linewidth. The magnetic permeability /i^ of CuCl is unity. Arrays of CuCl 
NPs can be fabricated by colloidal crystallization |4DJ and ion implantation techniques 
|41j . The CuCl NPs are coated with a metal of nanometre thickness, i = 0.105 = 
2.8 nm (nanoshell). Such hybrid metal-semiconductor nanoparticles have already been 
synthesized in the laboratory (for a recent review see Ref. |12]). The dielectric function 
of the metal is assumed to be described by the Drude model, i.e., 

e^ = l- , " (10) 
uj[uj + 17) 

In order to achieve NRI, ujp is taken to be Up = 1.05u;o. For the loss factor 7 we have 
taken a typical value of 7/0;^ = 0.01 [8J. 

The effective permittivity eg// and permeabihty /ie// of such a structure can be 
calculated by the extended Maxwell- Garnett theory [33], HI] which encompasses elements 
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of the Mie theory within the formulae of the eg//, /ie//- As such, the extended Maxwell- 
Garnett theory agrees very well with more rigorous approaches [HI [IHl [HI [IH [IS]- 
Figure [2^ shows the real and imaginary parts of Ce//, /ie//, and Ueff for the above 
described system of metal-coated CuCl nanoparticles. One clearly observes a region 
of negative 3ftee// for cu/cuo < 1-0015. One can also identify a narrow frequency region 
around w/ujo ~ 0.9998, where ^/J^eff < which is entirely within the region of negative 
^€eff- Within this region, Crieff < and a NRI band occurs. Note that for the 
calculation of the Ue/f the imaginary parts of eg//, /ie// have been also taken into 
account by choosing the value of n^ff = i/ee/zAie// which possesses positive imaginary 
part (the structures under study are made from passive materials). 

In order to verify the validity of the effective medium parameters depicted in 
figure [2K,b we have also employed the layer-multiple scattering method [3ll [35], [36] 
(briefly outlined in section [2]) in order to calculate the transmittance of light incident 
normally on a finite slab of an fee crystal of the above nanospheres. We have chosen the 
fee lattice since this type of lattice can satisfy a close-packed arrangement of spheres. 
The slab consists of 4 (001) fee planes of nanospheres and the respective transmittance 
is shown in figure [2b. It is evident that within the NRI of figure [2b, the transmittance 
shows a maximum which is what one expects from a slab of NRI [3]. As light propagates 
through the slab, it is also attenuated due to the intrinsic losses of the constituent 
materials (CuCl and Drude-metal) as it is evident from the corresponding maximum of 
the absorbance curve of figure [2t within the NRI band. It is also evident that there is also 
a wider transmittance peak from cu/cuq ~ 1.0015 to 1.002 which, however, corresponds to 
a region of refractive index with positive real part (the matching of these two frequency 
regions is not perfect as figure [2^ is based on an effective medium approximation whilst 
figure [2]d to an exact theory). In between these peaks as well as above the second peak, 
the transmittance is suppressed as only one of 3fJee// and 3^/ie// becomes negative giving 
rise to a practically imaginary refractive index (see figures [2bL,b). 

In order to probe the imaging properties of the above systems, we have applied the 
formalism developed in section [2] for the case of a single plane of metal-coated CuCl 
nanospheres. As a localized source of light we have considered a linear, infinitely thin, 
center-fed antenna of width d = O.OOIc/cuq, lying along the 2;-axis. It is placed at a 
distance h = 0.6c/ uq from the center of a single plane of spheres (see figure [I]). Figure [3] 
shows the distribution of all three electric-field components along the line y = in the 
image plane (see figure [I]), for frequency uj/ujq = 0.99984 which corresponds to the low- 
frequency maximum of the transmittance curve of figure ^jp. For comparison, we also 
show the corresponding curves in the absence of the plane of spheres. We note that the 
field distributions of all figures that follow contain both the far-field [g^ > (ky + g)^ in 
equation 1^] and near-field [g^ < (ky+g)^ in equation 1^] components. Since the chosen 
frequency {ujujo = 0.99984) lies within the NRI band of the array of spheres (see above) 
we expect that the near-field components which are transmitted through the plane will 
be amplified. And this is indeed what we observe in figure [3l We observe that all three 
components components of the electric field are amplified with respect to the case where 
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the plane of spheres is absent. This can occur only in the case of NRI where, as predicted 
by Pendry [3] , both polarization modes are amplified when transmitted through a planar 
NRI slab. We remember that when only one of eg// < or /Xg// < occurs, the near- 
field amphfication affects one of the polarization modes [2Sl ESI [2H1 ISH]- In addition 
to the near-field amplification of the electric field, we observe that the FWHM for the 
Ex— (~ Ao/7) and i^^— components (~ Ao/3) is considerably smaller than the free-space 
case suggesting a subwavelength imaging operation of the structure (see below). The 
E'z-component (figure [3b) is much more amplified than the other two; however, the 
field distribution looks much more complicated, possibly, due to interference effects (the 
far-field is also included in the calculation) and the peculiarity of the source (the Ez- 
component emitted by a center-fed antenna possesses a nodal line along the direction 
of the antenna axis). 

In order to confirm the fact that the near-field amplification shown in figure [3] 
is attributed to the occurrence of NRI, in figure H] we show the same quantities, 
for the same calculation setup as in figure [3l but for frequency u/uq = 1.00168 
corresponding to the high-frequency transmission peak of figure ^jp which lies out of 
the NRI band. We observe that only the i^^-component (figure Hb) shows a similar 
trend as that of figure [3b. The i^^-component is suppressed compared to the free-space 
case whilst the i^^.- component is barely amplified from the plane of nanospheres. Also, 
the subwavelength profile of the field distributions is lost as the maxima are almost as 
wide as those of the free-space case. So, judging from figure [H we can infer that the 
field amplification for all polarizations and the subwavelength extension of the electric 
field of figure [3 are in accordance with the occurrence of NRI. 

Next, in figure [5], we probe the imaging properties of thicker slabs, i.e., slabs 
containing more than one planes of nanospheres. We have chosen to depict the E^- 
component but similar effects are observed for the other two components. For the case 
of two planes of nanospheres (figure [5bi), the field amplification is still evident although 
somewhat distorted at the center. There are also two satellite peaks evident which 
may infiuence a potential lensing application of the structure. For a 3-planes-thick slab 
(figure [5]d) the amplification of the transmitted field is not as dramatic as for thinner 
slabs (Figs. [HH and [5bi) whilst for the 4-planes-thick slab the amplification is evidently 
lost. At first glance, this seems to be an unexpected result: the extended Maxwell- 
Garnett theory [131 HI] which was employed for the calculation of n^jf in figure [lb- 
refers to 3D collections of scatterers and it is is supposed to be more applicable to thick 
enough slabs. However, as both constituent materials of the nanospheres (metal and 
semiconductor) suffer from losses (especially the metal) we expect that thicker slabs 
mean longer optical paths for the transmitted field, and, therefore, higher losses which 
counterbalance the near- field amplification effect. This is also manifested in figure [2b 
where almost 90% of the incident power is absorbed within the slab (of 4 planes of 
spheres). Similar effects have been also measured experimentally for NRI metamaterials 
in the microwave regime [2^ . 
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Finally, in order to study the image-resolution properties of the array of 
nanospheres, in figure [HI we have considered the case of two localized sources (small 
center- fed antennas) separated by a vertical suhwavelength distance of c/ujq = Aq/vt. 
The array consists of a single plane of nanospheres since absorption is expected to 
deteriorate the image resolution for thicker slabs. The sources are placed symmetrically 
with respect to the center of a sphere, at a distance h = 0.6c/ Uq from the plane. As 
it is evident from the dashed curve of figure El the distance between the sources is too 
small to be resolved in free-space (a single peak appears). However, in the presence of 
the plane of nanospheres, the occurrence of NRI is accompanied by a dramatic increase 
in the image resolution which allows for distinguishing the two sources in space by use 
of visible light. We note that the image resolution can be further improved by using 
smaller nanospheres, down to radius of 15 nm. However, for CuCl spheres smaller than 
that, Heff is no longer negative [15] around the exciton resonance and a NRI band is 
not expected to occur. 

4. Conclusion 

It has been shown that arrays of metal-coated semiconductor nanospheres can exhibit 
NRI around the exciton resonance (if such exists) of the semiconductor. The existence 
of NRI has been studied by use of the extended Maxwell- Garnett theory and confirmed 
by rigorous EM calculations of light transmittance. Within the NRI band, an array 
of such spheres can amplify the near-field for both polarization modes, in accordance 
with Pendry's theory [3]. The near-field amplification is accompanied by suhwavelength 
image resolution that is only limited by the inherent losses of the metal and the 
semiconductor. The systems studied in this work are strong candidates for NRI 
metamaterials in the optical regime due to their fabrication with existing methods of 
colloidal chemistry |iT| liU] . 

References 

[I] Lakhtakia A and Mackay T G 2007 Opt. Photonics News 18 32 
[2] Veselago V G 1968 Sov. Phys.-Usp. 10 509 

[3] Pcndry J B 2000 Phys. Rev. Lett. 85 3966 

[4] Negative-Refraction Metamaterials, edited by G. V. Eleftheriades and K. G. Balmain (Wiley, New 
Jersey, 2005). 

[5] O'Brien S and Pcndry J B 2002 J. Phys.: Condens. Matter 14 4035 

[6] HoUoway C L, Kuester E F, Baker- Jarvis J and P. Kabos 2003 IEEE Trans. Antennas Propag. 51 
2596 

[7] Huang K C, Povinelli M L and Joannopoulos J D 2004 Appl. Phys. Lett. 85 543 
[8] Yannopapas V and Moroz A 2005 J. Phys.: Condens. Matter 17 3717 
[9] Felbacq D and Bouchitte G 2005 New J. Phys. 7 159 

[10] Wheeler M S, Aitchison J S and Mojahedi M 2005 Phys. Rev. B 72 193103 

[II] Wheeler M S, Aitchison J S and Mojahedi M 2005 Phys. Rev. B 73 045105 

[12] Jylha L, Kolmakov I, Maslovski S and Tretyakov S 2006 J. Appl. Phys. 99 043102 
[13] MacKay T G and Lakhtakia A 2006 J. Appl. Phys. 100 063533 



Subwavelength imaging of light by arrays of metal- coated semiconductor nanoparticles...9 

[14] Yannopapas V and Vitanov N V 2006 Phys. Rev. B 74 193304 

[15] Yannopapas V 2007 Phys. Rev. B 75 035112 

[16] Rockstuhl C, Lcdcrcr F, Etrich C, Pertsch T and Scharf T 2007 Phys. Rev. Lett. 99 017401 
[17] Peng L, Ran L, Chen H, Zhang H, Kong J A and Grzegorczyk T M 2007 Phys. Rev. Lett. 98 
157403 

[18] SchuUer J A, Zia R, Taubner T and Brongersma M L 2007 Phys. Rev. Lett. 99 107401 

[19] Redl F X, Cho K- S, Murray C B and O'Brien S 2003 Nature (London) 423 968; Saunders A E and 

Korgel B A 2005 ChemPhysChem 6 61; Shevchenko E V, Talapin D V, Kotov N A, O'Brien S and 

Murray C B 2006 Nature (London) 439 55; Kalsin A M, Fialkowski M, Paszewski M, Smoukov S K, 

Bishop K J M and Grzybowski B A 2006 Science 312 420 
[20] ParazzoH C G, Greegor R B, Li K, Koltenbah B E C and TanieUan M 2003 Phys. Rev. Lett. 90 

107401 

[21] Houck A A, Brock J B and Chuang I L 2003 Phys. Rev. Lett. 90 137401 
[22] Lagarkov A N and Kissel V N 2004 Phys. Rev. Lett. 92 077401 
[23] Grbic A and Eleftheriades G 2004 Phys. Rev. Lett. 92 117403 

[24] Aydin K, Buhi I and Ozbay E 2007 Appl. Phys. Lett. 90 254102; Ozbay E and Aydin K 2007 
Photon. Nanostruct.- Fundam,. Applic., doi:10.1016/j. photonics. 2007. 09. 003 

[25] Wiltshire M C K, Pendry J B, Young 1 R, Larkman D J, Gilderdale D J and Hajnal J V 2001 
Science 291 849 

[26] Wiltshire M C K, Pendry J B, and Hajnal J V 2006 J. Phys.: Condens. Matter 18 L315 
[27] Ono A, Kate J and Kawata S 2005 Phys. Rev. Lett. 95 267407 

[28] Fedorov G, Maslovski S I, Dorofeenko A V, Vinogradov A P, Ryzhikov I A and Tretyakov S A 

2006 Phys. Rev. B 73 035409 
[29] Alitalo P, Simovski C, Viitancn A and Tretyakov S 2006 Phys. Rev. B 74 235425 
[30] Yannopapas V 2007 Phys. Stat. Sol. (RRL) 1 208 
[31] Jackson J D 1975 Classical Electrodynamics (Wiley: New York) 
[32] Naus H W L 2007 Phys. Rev. E 76 026602 

[33] Modinos A, Yannopapas V and Stefanou N 2000 Phys. Rev. B 61 8099 
[34] Stefanou N, Karathanos V and Modinos A 1992 J. Phys.: Condens. Matter 4 7389 
[35] Stefanou N, Yannopapas V and Modinos A 1998 Comput. Phys. Commun. 113 49 
[36] Stefanou N, Yannopapas V and Modinos A 2000 Comput. Phys. Commun. 132 189 
[37] Gantzounis G and Stefanou N Phys. Rev. B 73 035115 

[38] Abramowitz M and Stegun I A 1965 Handbook of Mathematical Functions (Dover: New York) 
[39] Artoni M, La Rocca G and Bassani F 2005 Phys. Rev. E 72 046604 
[40] Orel Z C, Matijevic E and Goia D V 2003 Colloid. Polym. Sci. 281 754 

[41] Fukumi K, Chayahara A, Kageyama H, Kadono K, Akai T, Kitamura N, Mizoguchi H, Horino Y, 

Makihara M, Fujii K, Hayakawa J 1999 J. Non-Cryst. Sol. 259 93 
[42] Cozzoh P D, Pellegrino T and Manna L 2006 Chem. Soc. Rev. 35 1195 
[43] Doyle W T 1989 Phys. Rev. B 39 9852 
[44] Ruppin R 2000 Opt. Commun. 182 273 
[45] Yannopapas V 2007 Appl. Phys. A 87 259 

[46] Gaponenko S V, Germanenko I N, Gribkovskii V P, Zimin L G, Lebed V Y and Malinovskii I E 
1993 Physica B 185 588 



Suhwavelength imaging of light by arrays of metal- coated semiconductor nanoparticles. . .10 

t h h i 




Figure 1. (Color online) Calculation setup. 
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Figure 2. (Color online) The imaginary (a) and real (b) parts of the effective 
parameters ?^e// (solid line), eg// (dashed line) and fJ-eff (dotted line) for a 3D array 
of close-packed metal-coated CuCl spheres of radius S = 28 nm and coating thickness 
£ = 2.8 nm. (b) Transmittance, reflectance and absorbance of light incident normally 
on a slab consisting of 4 (001) fee planes of the above nanospheres, as calculated by 
the layer- multiple scattering method [311 ISSl [3S] ■ 
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Figure 3. (Color online) Distribution of (a), \Ey\'^ (b), and \E^\'^ (c) (in 

arbitrary units) along the line y = in the image plane (see figure [l]) , for frequency 
uj/oja ~ 0.99984. The solid lines refer to the field distribution from a localized source 
(antenna) placed at a distance h = O.Gc/wq from the center of a single plane of close- 
packed metal-coated CuCl nanospheres (see figure [T]). The dashed fines refer to the 
free-space case (no spheres present). 




Figure 4. (Color online) Same same as in figure[3]but for frequency oj/uiq = 1.00168. 
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Figure 5. (Color online) Solid lines: distribution of \Ex\'^ (in arbitrary units) along 
the line y = in the image plane (see figure [TJ, for frequency uj/luq = 0.99984, in the 
presence of a slab consisting of 2 (a), 3 (b) and 4 (c) planes of close-packed metal-coated 
CuCl nanospheres. The distance of the source from the center of left plane of the slab 
is h ~ 0.6c/ujo. The dashed lines refer to the free-space case (no spheres present). 




Figure 6. (Color online) Field distribution for two sources separated by distance 
c/luq = Xq/tt and arranged symmetrically with respect to the center of a sphere of a 
single plane of close-packed metal-coated CuCl nanospheres (solid lines). The sources 
are placed at a distance h = 0.6c/wo from the center of a single sphere of the plane. 
The dashed lines refer to the free-space case (no spheres present). 



